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Gravitational microlensing of quasar Broad Line Regions: 
the influence of fractal structures 
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ABSTRACT 

Recent models for the emission clouds within the Broad Line Region of quasars sug- 
gest that they are due to transient overdensities within an overall turbulent medium. 
If this were the case, the broad line emission would spatially appear fractal, possessing 
structure on a range of scales. This paper examines the influence of such fractal struc- 
ture when a quasar is microlensed by a population of intervening masses. It is found 
that while the highest fractal levels can undergo significant microlensing magnifica- 
tion, when these light curves are superimposed to create an emission line profile, the 
resultant emission line profile remains relatively constant for physical models of the 
Broad Line Region. It is concluded that the detection of the possible fractal structure 
of Broad Line Regions via gravitational microlensing is not practical. 
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1 INTRODUCTION 

While quasars are amongst the most luminous objects in the 
Universe, the majority of their emission is generated within 
a region only parsecs in extent. At cosmological scales, such 
regions are well below the resolving power of even the most 
powerful telescopes, and so the spatial structure within these 
inner regions cannot be directly observed and must be in- 
ferred by other, more indirect means. 

In recent years, a general picture for quasar central re- 
gions has emerged, with a hot accretion disk surrounding a 
supermassive black hole being responsible for the extensive 
continuum emission characteristic of quasars. Beyond this, 
on the scale of ~ lpc, this accretion disk is surrounded by 
a large population of clouds which are illuminated by the 
central engine and produce the broad emission lines, with 
widths of thousands of kms~ , also characteristic of quasar 
spectra 1 . 

Initial studies of the Broad Line Region (BLR) 
were based on simple ionisation models, with predicted 
BLR sizes of 0.1 to a few parsecs in extent (e.g. 
iDavidson fc Netzedll979h . Recent studies, utilising rever- 
beration mapping to measure the scale of the BLR, how- 



It should be noted that alternative models for the source of 
the broad line emission exist, including the accretion disk wind 
model of iMurrav et ail Il995ft which occurs on a considerably 
smaller scale. 



ever, have found evidence for a smaller BLR, more than 
an order of magnitude smaller than suggested by ionisa- 
tion models l| Peterson 1 Crenshaw 1 fc Meyersll98 5|). Further- 
more, IWandel. PeTeTsoTi^^Maikanri^99l ) and lKaspi et alJ 
( 2000) have demonstrated that the size of the BLR scales 
with quasar luminosity, with Rblr cx L ' 7 , and BLRs were 
also found to possess significant ionisation stratification, 
with high-ionisation lines arising in a region an order of mag- 
nitude smaller than low-ionisation lines. 

The origin of the BLR has proved problematic to 
understand, with suggestions that the emitting mate- 
rial could be the deb ris from star-disk collisions at 
the heart of the quasar flZurek. Stemk qnowska. fc Colgate! 
Il994l lArmitaee. Zurek. fc Davie s|^9^ , supernova exph> 
sions within quasar outflows ~ lPrttard et alJ l200ll 12003ft . 

or even the bloate d atmospheres of irradiated stars 

iAlexander fc Netzerlll997l) . However, the BLR must pos- 
sess significant substructure, with line emission arising 
from dense clouds embedded in a lower density medium 
JCapriotti. Foltz. fc Bvardl Il98lfl . with the smoothness of 
emission lines s uggesting that th ere must be of order 10 s 
emitting clouds iArav et al.lll998T) . effectively ruling out the 
picture of the BLR being composed of a relatively small 
number of large clouds or bloated stars. Nevertheless, such 
a situation is physically problematic, as these clouds should 
rapidly dissipate, with various mech anism s sugg ested for 
their confinement [e.g. magnetic fields llteesl il98it) l. 

Recently. iBottorff fc Ferlandl (1200 ll) suggested that the 
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'clouds' within the BLR are not isolated, individual enti- 
ties embedded within a confining medium, but rather rep- 
resented transient knots of higher density within an overall 
turbulent BLR. One conclusion of this study was that the 
BLR would not appear as a smoothly-varying emitting re- 
gion. Instead, the turbulent nature of the BLR would result 
in the spatial emission from the region possessing an overall 
fractal structure. 

Gravitational microlensing provides an opportunity to 
'see' the structure in the central regions of quasars, utilising 
differ ential mag nification effe cts llWambsganss fc Paczvnskil 
119911) . Recently. iLewis! (I2004D demonstrated that small-scale 
fractal structures within quasars, arising in an x-ray emit- 
ting hot corona above the accretion disk, imprint themselves 
on the light curve of a microlensed quasar. Hence, this paper 
considers the influence of microlensing on an extended, frac- 
tal BLR. However, as this is more extensive than the X-ray 
emitting region, it is important to compare it to the natu- 
ral scale length for gravitational microlensing; the Einstein 
radius. In the source plane, this is given by 



GM DuD oa 



D i 



(1) 



where M is the mass of the microlensing body, and Dy are 
angular diameter distances between the observer (o), lens 
(I) and source (s). For typical cosmological lensing configu- 
rations, with microlensing stars with mass M ^ 1M , this 
scale length is r\ ,$ O.lpc, and is often substantiall y smal ler. 
The X-ray emitting region considered in ILewis! l)2004h is 
small compared to this scale and hence signifi cant microlens- 
ing ma gnification is expected. As discussed in lLewis fc Ibatal 
(2004), however, the extensive nature of the BLR ensures it 
must lie across a substantial portion of the complex caustic 
network that is seen at high optical depths, and the effects of 
microlensing cannot be approximated by the influence of in- 
dividual caustic structures. The structure of this paper is as 
follows; Section[3]describes the numerical approach adopted 
in this study, while Section |3] discusses the resultant mi- 
crolensing light curves and their statistical properties. The 
conclusions to this study are presented in Section 2] 



2 METHOD 

2.1 Cloud Distributions 

The argument that BLRs possess fract al structure was ex- 
pound ed by the turbulence model of IBottorff fc Ferlandl 
(2001), and it is this procedure for defining a fractal cloud 
distribution that is utilised in this paper. The fractal struc- 
ture occupies a region with an overall scale size of Rmax- 
Within this radius there are a series of hierarchies whose 
radius is given by 



R(h) = Rr, 



(2) 



where h = 0, 1 • ■ • H is the level of the hierarchy, H denotes 
the maximum fractal hierarchy in the structure and L is a 
geometric factor. Each hierarchy possesses N substructures 
(the multiplicity) and the overall fractal dimension is defined 
to be 
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Figu re 1. A rea l isatio n of a fractal BLR as described in 
IBottorff fc Ferlandl 1200 J) . Clearly clumps and subclumps can be 
seen, but the smallest structures in this picture (~ 10 — 6 pc) are 
too small to be seen by eye. 



To generate a fractal distribution, a region Rmax is defined; 
this is the zeroth level of the hierarchy with h — 0. Within 
this region, N locations are chosen at random and represent 
the centres of the first level of the hierarchy (with h = 1) and 
assigned a radius given by Equation Within each individ- 
ual subregion, a further N centres are scattered at random 
and are assigned a radius corresponding to the second level 
of the hierarchy, with h — 2. The process is continued until 
the highest level of the hierarchy, h = H is reached; each 
point in this upper-most hierarchy corresponds to the lo- 
cation of a cloud with a radius given by Equation This 
procedure results in a total of n to t = N H clouds being dis- 
tributed over a region. For the purposes of this study, each of 
these final clouds was assumed to have the same luminoisty, 
although a more realistic model would probably require a 
distribution of luminosities throughout the region. 

In considerin g a realistic model for the BLR 
IBottorff fc Ferlandl (l200ll) determined the fractal param- 
eters required to reproduce known physical properties of 
the region, such as its column density and covering fac- 
tor. These values are employed in this current study, with 
Rmax = 5.12Rblr, where Rblr is the nominal radius to the 
BLR, L — 3.2 and N = 14 (see footnote 2). The maximum 
hierarchy which will be employed is related to the resolution 
of the microlensing magnification map, which is discussed in 
Section TFZ\ 

In choosing Rblr, the scale relation between 
the BLR radius and the luminosity of the quasar, 
as derived from reverberation mapping is empl oyed 
dWandel. Peterson, fc Malkanlll999l: iKaspi et all 120001) . As 
discussed in the next section, this study will focus on the 



D 



log N 
logL 



(3) 



1 In the study of lBottorff fc Ferlandl 120011) . N = 14.62, and frac- 
tal structures could be generated by drawing the number in a par- 
ticular hierarchy with this being the mean value. Fixing TV = 14 
does not significantly influence the resulting fractal structure. 
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Einstein radii 

Figure 2. The magnification map employed in this study. Rep- 
resenting image A of the microlensed quasar, Q2237+0305, a di- 
mensionless surface density of k = 0.36 and external shear of 
7 = 0.41 were employed. The units are in Einstein radii for a 
solar mass star and the scale of the image matches that presented 
in Figure HI 

microlensed images of the quasar Q2237+0305, which pos- 
sess an absolute magni tude of M ~ — 26 [accounting for a 
magnification of ~ 16 jSchmidt^Veb ster. fc Lewis! 1 1998 )]; 
as shown in llLewis fe Ibatal l2004) this corresponds to a 
Rblr ~ 0.05pc for high- ionization emission (e.g. C IV ) 
and Rblr ~ 0.4pc for low-ionisation emission (e.g. Mg II 
). For the purposes of this study, Rblr = O.lpc, a value 
between the high and low emission scale sizes. Using these 
parameters, a typical fractal BLR can be generated, with an 
example presented in Figure 

2.2 Gravitational Microlensing 

As noted in Section^ the substantial size of the BLR implies 
that the influence of gravitational microlensing cannot be 
approximated as a simple point-mass lens or isolated caustic 
structure. To this en d, t he backward ray-tracing ap proach of 
iKavser et al.l jl986t) and lWambseanss et al.l dl99fJ) was used 
to generate the magnification map employed in this study. 

The macrolensing parameters were chosen to represent 
Image A of Q2237+0305, the most intensively-studied mi- 
crolensed quasar, with a dimensionless surface mass density 
of k — 0.36 and an external shear of 7 = 0.41; this surface 
mass density describes the focusing due to matter within the 
beam of light traversing the galaxy, whereas the shear is the 
distorti on of the beam due t o the larger scale distribution of 
matter iKavser et al.lll986l) . All of the mass is distributed 
in microlenses with masses of 1 Mq . Given that the lens- 
ing galaxy in this system possesses a redshift of z; = 0.0395, 
while the source quasar is at z 3 — 1.695, the Einstein Radius 
(Equation for a solar mass star is r] ~ 0.06pc, adopting 
the current concordance cosmology. To accommodate the 
model BLR described in Section 12.11 the entire ray trac- 



ing region is taken to be 28 Einstein radii (~ 1.68pc) on a 
side. The resulting magnification map used in this study is 
presented in Figure 

2.3 Image pixelation 

Clearly, the fractal clouds outlined in Section 12.11 possess 
structure on a range of scales and this must be reflected 
in the magnification map. Hence the magnification map 
utilised in this study possesses a pixel scale of 1000 pixels 
per Einstein Radius, with one pixel physically correspond- 
ing to 6 x 10 _5 pc, and the total map presented in Figure [5] 
is 28000 pixels on a side. If a pixel in this map represents 
a single cloud in the fractal BLR, then the maximum hier- 
archy in the fractal distribution corresponds to H ~ 7 and 
the total number of clou ds in the BLR is Ntot ~ 1 .5 x 10 9 . 
Note that the model of iBottorff fc Ferlandl feOOlh consid- 
ered a maximum fractal hierarchy of H = 11, with a total 
of Nrot ~ 4 x 10 14 individual clouds; with this, the scale 
of each cloud would be ~ 1.4 x 10 _6 pc and the magnifica- 
tion map presented in Figure [5] would have to be 1.2 x 10 6 
pixels along each side, a very challenging prospect. As dis- 
cussed later, the lower resolution map adopted here does not 
significantly influence the conclusions of this paper. 

2.4 Generating light curves 

In typical gravitational microlensing problems, magnifica- 
tion maps are convolved with the surface brightness profile 
of a preferred source, providing a map from which the mi- 
crolensing light curve from that source can be drawn. Such 
a convolution, utilising Fourier transform techniques, was 
found to be computationally too expensive, given the num- 
ber of pixels used in the magnification map. Hence, in this 
present study, light curves were simply generated by extract- 
ing a single column of pixels along a magnification map at 
the location of a BLR cloud. Summing each of these columns 
for a collection of clouds together then gives the overall light 
curve of the BLR region of interest. With this approach, 
each individual BLR cloud (the smallest structure within 
the entire BLR) has a size of a single pixel in the map, cor- 
responding to 6 x 10 _5 pc. 



3 RESULTS 

Figure |3] presents a series of light curves derived from this 
analysis. Each panel presents subsections of the fractal hi- 
erarchy, from h = 7 in the top-left, consisting of a total of 
14 sources), to h — in the bottom-right, with a total of 
1.5 x 10 9 clouds. The background of each panel presents a 
greyscale map of the fractal structure (note in h — 7 and 
h — 6 the structure is too small to see). Each subpanel 
presents the microlensing light curve for the fractal hier- 
archy; the x-axis is in units of Einstein radii, whereas the 
y-axis corresponds to the magnification of the region. The 
expected crossing time of an Einstein radius in this system 
(assuming a tra nsverse velocity of th e lens of ~ 600km/s) 
is ~ 8 yrs (e.g. iLewis fc rrwinlli996T) . The degree of vari- 
ability changes quite strongly with fractal hierarchy h — 7 
displaying the strong variability, with caustic crossings quite 
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Figure 3. Example light curves for the fractal hierarchies presented in this paper, from the smallest (h = 7 top-left) to the largest (h = 
bottom-right). Each panel presents the cloud distribution (note that the smallest structures are not visible on the scale of these figures), 
with the resulting light curve arising from the magnification map presented in Figure [5] Each light curve is 1 Einstein radii in extent, 
with the location of the clouds in the figure representing the centre of the position at the centre of the light curve. 
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Figure 4. The magnification distributions for 250 realisations of 
fractal BLRs combined with the magnification map displayed in 
Figure |21 

similar to point-like sources 2 . Clearly these source sizes at 
fractal hierarchies of h > 7 are substantially smaller than 
the fractal structure seen in Figure [5] and hence the result- 
ing light curves for individual sources at this level would 
appear qualitatively the same as the h = 7 case, although 
the peak magnifications would be larger. Considering lower 
values of h, it can be seen that the light curve fluctua- 
tions ar e smoothed out, as expected as the source size in- 
creases iWambsganss fc Paczvnskilll99ll) . while at the base 
fractal level (h = 0), which considers the entire cloud popu- 
lation, the light curve is quite flat. 

To examine this further, 250 light curve samples were 
generated with differing realisations of the fractal hierar- 
chy. Figure 0] presents the magnification probability distri- 
butions for each of the fractal hierarchies in the entire sam- 
ple. These confirm the properties of the light curve seen in 
Figure |3] with the highest fractal hierarchy displaying mag- 
nifications of 10 — 20%. This distribution is quite different 
to the lower fractal hierarchies which display typical light 
curves for small sources. Here the magnification probabil- 
ity distribution is dominated by a peak at a magnification 
of 0.5, corresponding to extensive periods of demagnifica- 
tion in the light curve, with a high mag nification tail repre- 
senting the periods of caustic crossing dWambseanssI Il992l 
iLewis fc Irwiniri995f) . 

A simple interpretation of this result, therefore, is that 
if BLRs do indeed possess the fractal structure outlined, 
then, on the whole, they should not be significantly magni- 
fied by gravitational microlensing. While this is true, what 
it actually implies is that the total flux in the emission line 
will not vary as the quasar is microlensed, but it must be 
remembered that as well as spatial structure, BLRs also pos- 
sess kinematic structure which is responsible for the width 

2 The slight fuzz visible in the light curve is due to low-level 
numerical noise in the magnification maps. Its presence does not 
influence the results of this study. 



of t he emission line [see earlie r stud ies bv iNemirofj jl988f) 

and lSchneider fc Wambsgansll Jl990M 

As noted in iBottorff fc Ferlandl fcOOllh however, the 
fractal BLR model they present is purely geometrical and 
does not consider the kinematic aspects of the BLR. For the 
purposes of this paper, a toy model for the kinematic proper- 
ties is considered, with clouds of a particular fractal hierar- 
chy assigned a velocity drawn from a Keplerian distribution 
with random projections (i.e. all clouds at a particular frac- 
tal hierarchy and above are assigned the same projected ve- 
locity). This velocity structured BLR was then microlensed 
and the light curves summed into kinematic bins. The re- 
sults of this exercise are presented in Figure |S] the left-hand 
panels present the emission line profile, comprised of 1000 
velocity elements, whereas the right-hand panel presents the 
light curves for sampled values of the velocity elements over 
one Einstein radius. The spread in the emission line profile 
represents the superposition of all emission lines along the 
one Einstein radius. The pairs of panels arranged vertically 
represent the the various velocity models; i.e. in the h — 7 
model, each of the 14 clouds within the highest fractal hi- 
erarchy are assigned the same Keplerian velocity, for h — 6, 
14 2 within the second highest are assigned the same Kep- 
lerian velocity, etc. Clearly, the h = 7 emission line profile 
displays very little variability. Even at this high resolution, 
each velocity element must be the sum of of the light curves 
of many small sources, averaging out to give an overall flat 
light curve. In moving to lower fractal hierarchies, the emis- 
sion line profiles display more and more variability. However, 
it is clear that bulk motions of the lowest fractal hierarchies 
result in very "choppy" structure, with the h — 2 emission 
line breaking down into a series of spikes. The smoothness 
of observed emi ssion lines sugges t that such a situation is 
unphysical (e.g. lArav et al]ll998l) and hence the kinematic 
properties must be defined by the higher fractal hierarchies. 
It can be concluded, therefore, that the fractal structure will 
be virtually impossible to detect by studying the emission 
line profiles of microlensed quasars. 



4 CONCLUSIONS 

The nature of the broad line emission regions of quasars 
has remained a subject of discussion for a number of years, 
with recent models suggesting that the emitting clouds are 
temporary density enhancements in an overall turbulent 
medium. Such a scenario predicts that these density en- 
hancements should have a fractal distribution of structures. 
This paper has investigated the influence of gravitational 
microlensing on this extensive fractal structure, finding that 
different levels of the fractal structure undergo differing mag- 
nifications, with the smaller sections of the substructure suf- 
fering stronger magnification. It was found, however, that 
while individual clouds in the fractal hierarchy displayed 
quite dramatic variability due to microlensing, the combi- 
nation of light curves from the entire population of clouds 
within the BLR resulted in an overall constant light curve. 

While the overall magnification of the BLR appears to 
be constant, more significant magnification of substructures 
could be apparent when examining the form of the emis- 
sion line profile, i.e. substructures within the BLR could 
possess coherent velocities and hence may contribute to a 
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Figure 5. The emission line profile for the toy kinematic models 
presented in this paper. The left-hand panel presents the emis- 
sion line profile for the appropriate bulk motion of the fractal 
structure. The right-hand panel presents the light curve for vari- 
ous portions of the lines, colour-coded to match the wavelengths 
presented in the left-hand panels. In this figure, the total fractal 
hierarchy is the same for all frames and the h values denote the 
hierarchy at which the bulk velocity has been assigned; i.e. in the 
top panel with h = 2, each h = 2 has been assigned a velocity 
drawn from the model and this value has been assigned to each 
higher fractal level within that h = 2 structure. Hence the h = 7 
possesses more kinematic diversity that h = 2, resulting in the 
overall smoothness of the resulting emission line. 



particular velocity within the emission line. This paper con- 
sidered a simple Keplerian velocity structure for the emis- 
sion line clouds, assigning a particular velocity to all fractal 
hierarchies below a particular level. The result of this pro- 
cedure, however, revealed that coherent velocity structure 
cannot apply to the lowest levels of the fractal hierarchy as 
the resultant emission line profile is clearly too structured, 
although dramatic variability is seen throughout the line. 
Assigning coherent velocity structure to the higher fractal 
hierarchy does smooth out the form of the emission line, 
but it also smooths out the light curves for each velocity 
bin. From this it is possible to conclude that while individ- 
ual structures within a BLR are being microlensed and are 
undergoing significant magnifications, the extensive nature 
of the BLR and the requirement that the resulting emission 
line appears relatively smooth means that microlensing is 
unlikely to reveal the putative fractal structure of quasar 
BLRs. 
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